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Abstract 

The notion of a Laplace ladder for a discrete analogue of the 
Laplace equation is presented. The adjoint of the discrete Moutard 
equation and a discrete counterpart of the nonlinear form of Goursat 
equation are introduced. 



1 Notation 

Value of functions of continuous variables we denote by f(u,v) i.e. / : M? 3 
(u, v) I— > f(u, while functions of discrete variables by /(mi, m 2 ) i.e. / : 

Z 2 3 (mi, 777.2) 1— > /(mi, 777.2) € K. Partial derivatives are denoted by comma 
e.g. f, uv (u,v) := g^(-u,f). Shift operators are denoted by subscripts in 
brackets e.g. f(m 1 ,m 2 ) {1) := /(mi + l,m 2 ), f{m u m 2 ) {2 ) ■= f{m 1 ,m 2 + 1), 
/(mi,m 2 )(i2) := /(mi + l,m 2 + 1), /(mi,m 2 )(-i) := /(mi - l,m 2 ) etc. We 
omit arguments when operators indicate what kind of functions (of discrete 
or continuous variables) we deal with. Difference operators are denoted by 
Aif := f\i) — /. The diamond operator we define as follows O/ := j^f ■ 

2 Introduction 

The Laplace equation 

Lijj(u,v) = (1) 
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where L is a differential operator 



L = & ' A(u, v)l- B(u, v)f v - C(u, v), (2) 

is often regarded as a master equation of theory of .S-integrable systems. 
One can justify this point of view indicating on either the existence of large 
class of Darboux-type transformations introduced by Jonas |I[ and Eisenhart 
(called fundamental transformation or Darboux-type transformation for 
multicomponent KP hierarchies) that acts on solution spaces of set of the 
Laplace equations and give rise Darboux-Backlund transformations for the 
Darboux equations or modern approach of 9-dressing method || that was 
applied to Darboux equations by Zakharov and Manakov. 

Both 9-dressing method and, what more important for us, discrete version 
of fundamental transformation acting on solution spaces of discrete counter- 
parts of the Laplace equations 

Cif)(m 1 ,m 2 ) = (3) 

C = A(i)A( 3 ) - A{m li m 2 )A ( i) - B(m 1 , m 2 )A (2) - C(m 1 ,m 2 ), (4) 

was applied to discrete analogues of Darboux equations 0, |[ ||. Moreover 
a theory of reductions of the discrete fundamental transformation is being 
developed 0, §, |, [l^, [11], [12 . 



On the other hand Athorne |13[ indicated that a symmetry of a Laplace 
ladder of Laplace equations impose constraints on the operators L of equa- 
tions of the ladder. In particular the cases when Laplace ladder contains a 
Moutard equation 

ip,uv= f(u,v)i> (5) 

and Goursat equation 

ip,uv= l(log\(u,v)), v iJj lu +\(u,v)4> (6) 
which is related to [T^, |1^ 

$ xy = 2y/\(x,y)0 x y (7) 
were discussed. There exist Darboux type transformations that leave the 



form of Moutard equation |T4| and Goursat equation |L5| [L6| unchanged. 

We start from the basic facts from the theory of discrete Laplace equation 
(the invariant description and T-equivalence are introduced in section [3D. 
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Next we introduce two operations acting on invariants of discrete Laplace 
equation: Laplace transformations (following the papers \ T7\ , [18 , p79f| we recall 
Laplace sequence of discrete Laplace equations ) in section ^ and an adjoint 
operation in section |5|. On combining the notions of Laplace sequence and 
adjoint operation we introduce in section || Laplace ladder of discrete Laplace 
equations. The Laplace ladder turns out to be a proper object to discuss an 
old and a new discrete Moutard type equation that respectively appear first 

(see section 0) and the discrete Goursat equation 



and 22, 



in the papers 
ill as well (see section 



3 Discrete Laplace equation. T-equivalence. 

A differential operator (P does not change its form under the transformation 

L^L = - g oLog (8) 

where g denotes an operator of multiplying by a scalar function g = g(u, v) ^ 
The above transformation is usually called gauge transformation. We say 
that two Laplace equations Ltp(u,v) = and Li/j(u,v) = are equivalent, 
iff' there exists a function g such that operators L and L are related by (§) 
(Relation: {two Laplace equation Li(j(u,v) = and Li/j(u,v) = are in the 
relation iff there exist a gauge transformation (Q)} is equivalence relation). 
The following functions (Laplace invariants) are invariant with respect to the 
gauge transformation 

h — AB — A, u +C, k = AB — B, v +C. (9) 

If the invariants of the operators L and L are equal (k = k and h = h) then 
there exist a gauge transformation that transforms L into L. Ordered pair of 
invariants (h, k) labels the equivalence classes of equations (|l|). We refer to 
Athorne paper [TBj for further information concerning the continuous case. 
The discrete Laplace equation @ can be written in the form 

^(12) - cr0(i) - Pi>(2) -j?p = Q (10) 

where 

a:=A + l (3:=B + 1 j:=C-A-B-l (11) 
By analogy to the continuous case we introduce invariants 

k:=^, n:=^. (12) 

7(2) ' 7(1) v ' 
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with respect to a discrete gauge transformations £ \— > j^~£<?> which is the 
discrete counterpart of the gauge transformation (§) and where g denotes 
operator of multiplication by function g = g(mi, 117,2) 7^ 0. We shall call the 
invariants basic invariants. They are basic in the sense that if two operators 
£ and £ have equal basic invariants ((n, k) = (n, k)) then they can be related 
by a gauge. 

For our convenience (section |7|) we introduce also derivative invariants 
which do not poses the basic property. Secondary invariants of the first kind 
are defined by 

/C := KK(i) 7i := nn^ 2 ) (13) 
Secondary invariants of the first kind are defined by 

K := ^2) H := ^ (14) 

ft n V / 

We denote equivalence classes of the discrete Laplace equations by (n, k) . 

If the basic invariants of two discrete Laplace operators £ and £ are 
related as follows 

k = Tk, n = Th, 

where T denotes a shift operator (i.e. Tf(mi, 777,2) = + k, ni2 + 1) where 
k G N and I G N) then we shall say that equations are T -equivalent and 
we denote the enlarged equivalence classes by T(n, k). T-equivalence plays 
important role in our considerations. 



4 Laplace transformations and Laplace sequence 

In order to solve a discrete Laplace equation (0) we can try to factorize the 
discrete Laplace operator (f|). There are only two ways we can do it 

(A (1) - B)(A (2) - Ai-dW - W = 0, (15) 
(A (2) - A)(A {1) - B ( _ 2) )vp - hj) = 0, (16) 

where 

(3 = C-A + A i - 1) (B + 1) =7 + , , 

! = C-/3 + % 2) (^+l) = 7 + A-2)«- 1 J 

If one of the functions f) or 6 is equal to zero we succeed in factorization (and 
we can solve equation by quadratures). If not we can introduce functions 

^ T := (A (2) - A { - X) )i) (18) 
^ := (A (1) - /3 ( _ 2) )^ (19) 
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equations flT5|) and flT6"| ) can be written in the form 

^ = i(A (1) ^ T -^ T ) (20) 

V; = |(A (a) ^-^) (21) 

If we eliminate from equations (|I8|-|T9|) the function ip using ([20j-f2T|) then we 
obtain that functions ^ and ip^ satisfy equations 

tPU = + + 7^ T (22) 

4=0 = a «4) + x^^) + ¥ 7 ^ (23) 

and invariants of these equations are connected with invariants of the equa- 
tion © via 



n - K V (l+n) ( _D ~ n ( 2 ) 



(24) 



(2) 



(25) 



4) = ^ = ^<>(^ 

H T = (HO(l + n)) H) K T W = (K (1) 0^ 

H ( 2 ) = K)°(lfe) (1) Ki = ( K °(! + *))(-*) 

One can easily show, that 

(n , /c ) = (n^ , k^) = (n, «;) 

so Laplace transformations for the discrete Laplace equations i.e. maps 
(n, k) i — ► and (n, k) h- > (rr, are mutually inverse. By analogy 

to continuous case we call the sequence of the equations 

(n , /c ), (?i , z^), (n, «), (n\ re"), (rr\ re^), ... (26) 

Laplace sequence of discrete Laplace equations. Three comments are in order. 
Firstly, Laplace transformation are well defined on the T-equivalence classes. 
Secondly equation 

J „l - 1 



n (l) n (2) = nn (12)<>J 



+ n 
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where the Laplace transformation j (J,) are treated as increment (decrement) 
in third discrete variable is nothing but a form of Hirota equation ||20|| . Finally 



if we introduce X := 4- and 3^ := - where a and 6 are a solutions of fe(i) = /36 



l (2) — «(_i)a then equations ( |T8D and take form 

A (2) ^ = 

A m x = vy 



(27) 



where 

= -*- (28) 

0(1) 0(2) 

and equations ( ^2]) and (pTOf ) in the new gauge are 

^(12) = y w + ^y { 2) - M^y (29) 

*d2) = ^X {1) + X (2) -M^-X (30) 

where 

M = 1-VQ. 



5 Fundamental transformation and adjoint 

Fundamental transformation that maps from solution space of equation ([I]) 
has two functional parameters @. The first one is a solution of equation ([I]) 
and the second one is a solution of adjoint of equation ([[]) 

^, uv + {Aft) , u + (B^) , v -Cy = 0. (31) 

Our aim is to introduce the discrete adjoint of equation (fL0|). To do it we 
shall suitable rewrite fundamental transformation P for the single discrete 
Laplace equation (JJ). 

Discrete fundamental transformation maps from the solution space of the 
equation (we choose for convenience the so called affine gauge i.e. C = 0) 

A (1) A (2) x = + -f-^2)x (32) 

to the solution space of the equation 

1 Afoia 1 , A(i)b l , 
A^A^x 1 = -^A^x 1 + -^A^x 1 (33) 
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and is given by 

The functional parameters of the transformation d i #' are not arbitrary, 
namely the function $ is a solution of equation ( |32D i.e. 



A (1) A (2) = ^A (1) i5> + ^A (2) tf (35) 
while function we construct from i? and a solution i?t of the equation 

A (I) A (2) # + A (1) ( 1+ A^ + i y j + A (2) {—^-^ j = (36) 

The way to obtain d' looks as follow: first we find auxiliary functions A i \ 

A wy A (2)° 1? t 

A ( i)A = - ■ A(2 ; a ■ A(i)6 A (2) A = -A (2) # _ A ^ a 

X = A + 0t 

and next function we are searching for 

A (1) tf' = X A (1) 0, 
A (2 )t?' = AA(2)0. 



(37) 



(38) 

Conditions (^-|3S|) assure us, that x 1 given by (133) satisfies (P3"[). New func- 
tions a 1 and 6 1 are related to old ones via 

A(2)»' _ A (2) (a X Jr-a) A ( i)b 1 _ A (1) (b\^-b) . , 

a 1 - (a X £-a) ' ~ (bX^-b) ^ J > 

We say that equation fl3"6| ) is adjoint of equation (^). Invariants n^, «t, 
H"!" and of equation ( |3~6l) are conected with invariant n, k, H and K of 
equation (|35|) in the following way 

n 1 " = «(x) «t = n(2)) (40) 

Ht = K ( i) Kt = H (2) . (41) 
We say that equation £^ = is adjoint of equation £^ = iff 

r(7it,«;t) = T(K (1) ,n (2) ). (42) 

We observe that T(n, = T(n, k). 
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6 Laplace ladder 

On superposing Laplace and adjoint transformations we get 



y-l jK ti) = I v— — .// 





K 


l+n ( 


-1) 


n(- 


l)n(2) < 








K 




-1) 


*(- 


2)«(1)< 


> 1 






n 




-2) 


*(- 


2)«(1) < 


> 1 





( n tT, K tT) = I - ^^ V T -f— , AC 



(43) 



(2) 
'(1) 
(1) 

s (-2) ' /(2) 

It means that we have 

T(ratt, /t Tt ) = T(nti, K ti), T(n^, K it) = 2>tt, K tT). (44) 
We obtain a commuting diagram Laplace ladder of discrete Laplace equations 

— ► T(n T ,K T ) — > T(n TT ,«; TT ) ... 

... T(ntT )K tT) < — r(r»t j/c t) < — T^,/^) < — T{n ni ,K Ul ) ■■■ 

where the long arrows denote f Laplace transformations while vertical lines 
(rungs of the ladder) the adjoint operation f. 

7 Moutard subcase 

A Laplace equation is called discrete Moutard equation iff its invariants are 
related by 

H = K, (45) 

or equivalently H = /C or n^n = ac(i)AC. Every discrete Moutard equation 
can be transformed to equation 

N (12) + N = F(N {1) + N {2) ), (46) 



by gauge transformation. The equation fl4"E| ) together with the Darboux type 
transformation has been introduced by Nimmo and Schief 0. 

From formulas (fy]) we get that an invariant characterization of the adjoint 
of discrete Moutard equation is 

R (2) ~ XV (1) 



H+ =Kt (47) 



So opposite to the continuous case the discrete Moutard equation is not 
self-adjoint. A Laplace equation is called adjoint discrete Moutard equation 
iff' its invariants are related by (f4Tp. Every adjoint discrete Moutard equation 
can be transformed by a gauge transformation to equation 

Af ( i2) + AT = F (1) JV (1) + F (2) M (2) . (48) 

If function iV satisfies a discrete Moutard equation ( pE6|) then function M 
given by 

A" := N {1) + N {2) (49) 

satisfies the adjoint discrete Moutard equation (||3|). Due to relationship ( j4Tf ) 
equation (f48"D "inherit" integrable features from equation (|46|) . In particular 
from the discrete Moutard transformation T7T1 we have 



(50) 



where 6 is solution of discrete Moutard equation ( [46] 

e (12) + = F(e (1 ) + e (2) ). 



Further features of the equation (f48|) will be discussed in p6 



8 Goursat (symmetric) subcase 

On demanding 

rJ = K (51) 

or equivalently 

H T = K (52) 

the equations (|T0D and (|2^), which are equivalent respectively to equations 
(l29| ) and (|30|), become the discrete Goursat equations. Invariant characteri- 
zation of the discrete Goursat equations is 



r (u lV ( 53 ) 

^(l)) 2 = ^(12)0^ 
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In terms of functions given in we have \'.\ I 



0: 



On introducing function r 



T (12) T 
r (l) T (2) 



M = l- VQ 



(54) 



(55) 



and due to the fact that the equations ( p7|) has the symmetry 
without lost of generality we can take 

VT(x) = Qt( 2 ). 

On defining 



x:=X 



we get pi| 



1/(2) 



(56) 
(57) 

(58) 



If functions x and y satisfy the above system then there exists a function 9 
such that 

x 2 = A (2) 2/ 2 = A (1) (59) 

We also have 

1-^ = ^21(1) = Q2IM 

T (2) r (l) 

and therefore squared equations (|58D (discrete analogue of (0)) take form 
1-M 



A(i)A(2)0 



(A,,,^ A l2 ^)±2\/i^ A /A (11 ^/A i; ^ (GO) 



where ± is result of square rooting of (|59|). The equation (|60|) we shall call 
a nonlinear version of discrete Goursat equation. 

In the case of Goursat reduction just as in the continuous case the ladder 
twist and fold to the diagram 

T(n\ri) — , . — , . — , . — , . ... 
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for in virtue of formulas (P%|), ([II]) and (^Tj) we have (n^, k^) = (/%), nJ 2 )) 
(w(2),AC(2)) and analogously (n^,K^) = i.e. 



T(n Tt , /« Tt ) = T(n, k), T{n\ k^) = T{n\ k t ). 

9 Conclusions 

We have shown that the Moutard equation splits in the dicrete case into dis- 
crete Moutard equation and adjoint discrete Moutard equation. This struc- 
ture should be taken into considerations, when a discretization of systems of 



a physical [24, 22, \Z3 and of geometrical [25, 26 1 importance is discussed. 
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